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Abstract 

A generalized Tomonaga-Schwinger equation, holding on the entire boundary of a 
finite spacetime region, has recently been considered as a tool for studying particle scat- 
tering amplitudes in background-independent quantum field theory. The equation has 
been derived using lattice techniques under assumptions on the existence of the con- 
tinuum limit. Here I show that in the context of continuous euclidean field theory the 
equation can be directly derived from the functional integral formalism, using a technique 
based on Hadamard's formula for the variation of the propagator. 

1 Introduction 

Quantum dynamics can be formulated in a variety of ways. One of these, championed by 
Feynman pQ, is to give the propagator of the Schrodinger evolution W(qi,q2, t). Here q\ 
and q2 are initial and final positions and t is the propagation time. The propagator satisfies 
the Schrodinger equation in the initial as well as in the final variables. Formally, it can be 
expressed as the path integral of the exponential of the action, over all paths that start at q± 
at time zero and end at q<i at time t 

W( qi ,q 2 ,t)= L n Dxe-rA-] . (l) 

J x(t)=q<2 

Quantum field theory too can be formulated in a variety of ways. One of these is the 
Schrodinger functional representation (for a review, see 0), where states are expressed as 
wave functionals of the field ip defined on a spacelike 3d plane. Evolution is given by a 
functional-differential Schrodinger equation. The propagator of this equation is a functional 
W[(fi,(f2, t], which propagates the initial field ipi at time zero to the final field <^2 at time 
t. I give the explicit form of W[<pi,(p2, t] for a free field theory in an appendix below. The 
functional W[(pi,(f2, t] satisfies the functional Schrodinger equation on both the initial and 
the final spacelike 3d planes. Formally, it can be expressed as the path integral of the expo- 
nential of the action, over all field configurations that take the initial value cpi on the initial 
plane and take the final value 922 on the final plane. 

Robert Oeckl 3. has recently argued that it should be possible to formulate quantum field 
theory in terms of a functional W[(p, S] of the field <f> defined on the compact 3d boundary £ 
of a finite spacetime region 7Z. In this case, £ is not going to be spacelike. The same point 
of view has been extensively developed by Carlo Rovelli in the book [3], to which I refer for 
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a detailed discussion. Formally, W[y>, E] can be expressed as a path integral over all field 
configurations on TZ, that take the boundary value <p on the 3d boundary E 



The field ip can be viewed as expressing initial, final as well as boundary values of <£, and 
W[(p, E] expresses the amplitude for having a certain set of initial and final fields, as well as 
boundary fields, measured by devices that are located in spacetime in the manner described 
by (the geometry of) the surface E. This picture is near to what actually happens in a 
laboratory experiment, where the initial and final state of a scattering event are confined into 
finite-size spacetime regions. The relation between particle states that can be defined in such 
a finite context and the usual particle states of quantum field theory, defined on an infinite 
spacelike region, is discussed in 

The interest of this suggestion is in its relevance for quantum gravity. Here, diffeomor- 
phism invariance (that is, background independence) implies that W[ip, E] is independent 
from E. The key observation, pointed out in [fy, is that the boundary value of the gravita- 
tional field on the boundary surface E determines the geometry of E, and therefore codes the 
relative distance of the measuring devices and the proper time elapsed from the bottom to 
the top of the spacetime region considered. Therefore in gravity the background independent 
propagator W[(f] contains the same physical information about the geometry of the bound- 
ary surface as the propagator W[<p, E] of a background dependent quantum field theory. In 
a scattering experiment we measure incoming and outgoing particles, that is, matter-fields 
variables, as well as distances between instruments and elapsed time, that is, gravitational- 
field variables. The two kind of variables are both represented by the boundary fields (p. In 
other words, the boundary formulation realizes very elegantly in the quantum context the 
complete identification between spacetime geometry and dynamical fields, which is the great 
Einstein's discovery. A simple model where this picture is concretely realized is presented 
in [7j. Thus, it was suggested in jS] that particle scattering amplitudes can be effectively 
extracted from W[<p] in a background independent theory; details will be given elsewhere (Hj. 

Given the possible relevance of the propagator W[<p, E] for quantum gravity, it is impor- 
tant to investigate it in detail, and to check whether its expected properties truly hold. In 
particular, it was argued in jS] that W[</?, E] should satisfy a local functional equation gov- 
erning the variation of W[<p, E] under arbitrary local deformations of E, namely an equation 
of the form 



Here s is a coordinate on E, and the precise meaning of the other symbols will be defined 
below. When E is formed by two parallel planes, © becomes the field propagator, and its 
variation under a parallel displacement of one of the planes is governed by the functional 
Schrodinger equation. The possibility of extending this equation to variations of arbitrary 
spacelike 3d surfaces was explored by Tomonaga and Schwinger already in the late forties 
jUj. These authors derived a well-known local generalization of the functional Schrodinger 
equation, called the Tomonaga-Schwinger equation. The functional equation (pj) generalizes 
the Tomonaga-Schwinger equation, since it holds for arbitrary boundary surfaces. Indeed, 
if the general philosophy of the boundary approach is correct, the distinction between the 




(2) 




(3) 



2 



initial and final fields on spacelike portions of E on the one hand, and the boundary values 
of the field on the timelike portions of E, should become of secondary relevance. 

Equation Q has been derived in ^0] on the basis of a lattice regularization of the func- 
tional integral © defining W[(p, E], and under certain hypotheses on the existence of the 
continuum limit. Here, working in the context of a free euclidean theory, I show that this 
equation can be derived from the functional integral definition of W[(p, E] directly in the 
continuum, using a formula by Hadamard which expresses the variation of a Green function 
under variation of the boundary 11 . Although incomplete, I think that this is a relevant 
step towards establishing the general viability of an equation of the form ©. 

Furthermore, as a preliminary step in the derivation of (J3J), I derive its classical limit, 
which is a generalized Hamilton-Jacobi functional equation of the form 



<5E(s) V M*) 

satisfied by the classical Hamilton function S[cp, E]. The Hamilton function S[(p, E] is the 
value of the action computed on the solution of the equations of motion that takes the value 
(p on E (see jl].) The generalized Hamilton- J acobi equation (jlj) is extensively discussed in 

a- 

The paper is organized as follows. In Section 2, I define the theory and equation (J3J); in 
Section 3 I introduce some Green function techniques; in Section 4 I prove equation ((IJ; in 
section 5 I prove equation ©. In the Appendices I sketch a proof of the Hadamard formula, 
and I give the explicit expression of the functional Schrodinger propagator W[tpi,ip2, t] for 
the euclidean and for the lorentzian theory. 



2 Definitions 

2.1 Surface and surface derivative 

Consider a finite region 1Z in the euclidean 4d space i? 4 . I use cartesian coordinates x,y,z ... 
on i? 4 , where x = (x' 1 ), fi = 1,2,3,4. Let E = dlZ^ be a compact 3d surface that bounds a 
finite region TZy,. I denote s,t,u. . . coordinates on E, where s = (s a ), a = 1,2, 3. Then E is 
given by the embedding 

E:s a ^a^(s). (5) 
The euclidean 4d metric, which raises and lowers the /U indices, induces the 3d metric 

on E, which can be used to raise and lower the three dimensional a indices. The surface 
gradient V is defined as 

V a = ^- (7) 

The normal one-form to the surface is 

dx"{s) dxP(s) dx°(s) 
n^s) - e^pv ggl gg2 gg3 . (»j 
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I orient the coordinate system s so that is outward directed. Its norm is easily seen to be 
equal to the determinant of the induced metric 



= det q. (9) 

In the following, I will use the normalized normal 

Ufj, = (detg)~^ (10) 

and the induced volume element on E 

dS(s) = (det<?(s))5 d 3 s. (11) 

(I shall simply write dT, when the integration variable is clear from the context.) Notice that 
the combination dT, does not depend on the metric, since 

d£(a) n^s) = d 3 s n^s). (12) 

Given a functional F[T] that depends on the surface, I define the functional derivative 
with respect to the surface as the normal projection of the functional derivative with respect 
to the embedding © that defines the surface 

n"(s)-=- ^ T -. (13) 



Here the functional derivative on the right hand side is defined in terms of the volume element 
c?E. That is 11 , 5F [S]/5S(s) is defined as the distribution that satisfies 

m Sf - m £$ - . m 

where S e jv is the surface defined by the embedding 

S £7V :s a i — >x^(s) +eN(s)n^{s). (15) 

Geometrically, this derivative gives the variation of the functional under an infinitesimal 
displacement of the surface in the normal direction. 

The surface derivative takes a simple form for certain simple functionals. Consider a 
scalar field f(x) on R 4 . Let Ff[E] be the integral of / in the region TZ^, bounded by S; that 
is 

F f [T] = [ d 4 xf(x) . (16) 

Then easily 

if - /«•»■ w 

That is, the variation of the bulk integral under a normal variation of the surface is the 
integrand in the variation point. Similarly, consider a vector field v^(x) on R . Let i^P] be 
the flux of v fJ, (x) across S; that is 

= J dE(a) n M («K(«) • (18) 
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Then Stokes theorem gives easily 

5F„m 

£E(s) 



d^{x{s)). (19) 



That is, the variation of the flux by a normal variation of the surface is the divergence of the 
vector field. 

2.2 Field theory and Tomonaga-Schwinger equation 

I consider a free euclidean scalar field 3>(x), defined in 7Z, with assigned boundary conditions 
<p(s) on E. That is 

= <p{s). (20) 
The dynamics of the field is governed by the lagrangian density 

C = + ^m 2 $ 2 , (21) 

and the equation of motion is 

(-□„. + m 2 ) $ (x) = 0, x e V (22) 

with the boundary conditions 1)20(1 on E. Here H\ x = d/dx^ d/dx 11 . By choosing one of the 
coordinates x M , say x 4 as a "time" coordinate, it is possible to derive the hamiltonian density 

F($,V<I>,n) = in 2 -i(V$) 2 -im 2 cI> 2 ; (23) 

where II is the canonical momentum associated to the field <J> and V is the gradient on the 
x 4 = constant surface. 

For this theory, the integral in is a well defined gaussian functional integral. Therefore 
W[(p, E] is well defined. I show in this paper that W [<p, E] satisfies the generalized Tomonaga- 
Schwinger equation 



5E(s) V M«) 

= (2 ("Mi)) -2^W) a -2 m V(.)J^,q- 

3 Green functions 

The tool I use to prove (|24j) is a Green- function technique. Given E, let the Green function 
Gs(x,y) be the solution of the inhomogeneous equation 

(-n x + m 2 )Gz(x,y) =5&(x-y) (25) 

in the region 1Z, satisfying the boundary condition 

G s (x(s),y) = , (26) 



5 



namely that vanishes when x is on E. I introduce the following useful notation: 

<9G s (x,y) 



x=x(s) 



(27) 



and similarly for d^Gj](y, s), d^d^G^is, t), and so on. The solution of (|22|). with the bound- 
ary condition 1)20(1 . namely with the boundary value tp on E, is then given by 



(28) 



%,x(y) = - / dE(s) <p(s) d s n Gx(s,y) ■ 
Jt, 

To see that this is indeed the solution, one can write the equality 

$ (y) = J <fx (* (x) [-U x + m 2 ) G E (x, y) - G s (x, y) {-U x + m 2 ) $ (x)) , (29) 

which is satisfied because of ()25(1 and (|22|> . The right hand side can easily be written as a 
surface integral, which gives 1)28(1 using the boundary conditions ()20|) and 1)26(1 . 

The propagator W[ip, S] can be written explicitly in terms of the Green function. In fact, 
it is possible to solve the gaussian integral (j2j) using the fact that the classical solution (f28|) 
is an extremal value of the exponent in ((2J). The integral is then a gaussian integral over the 
fluctuations. The action on the classical solution with given boundary conditions is called 
the Hamilton function |12j : 



S[ip,T] = S[%^] 



2 JVa 



(fx 



(30) 



Inserting 1)28(1 . it is possible to write this as 

1 



S[ip, S] 



(31) 



It is then easy to perform the gaussian integration in (J2J), which gives 
exp(-5[* v ,s]) 



W[<p,E] 



y/det (-□ + m 2 ) 



VdetGs (x,y) exp ( - y dE(«)dE(t) <p(a) <p(t) d^G^t) 



(32) 



In the following, it is necessary to have an expression for the variations of Green function 
Gy, with respect to a displacement of the surface S. This can be done using Hadamard's 
formula |llj 

5Gt, (x,y) 



<5£(s) 



d s n Gx(s,x) d s n Gx(s,y) . 



(33) 



I sketch a proof of this formula in the Appendix A. 
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4 Hamilton— Jacobi 

In this section I prove that the Hamilton function (|31j) satisfies the generalized Hamilton- 
Jacobi equation (£Q). From (|2*H|) . this reads 

+ 2 l^rJ - 2 m ^ (s) - 2 (v ^ (s)) = • (34) 

To calculate the first term of this expression, I use the form (|3U|) of the Hamilton function 
5S[ip,E] 8 If, 



( d 4 x Ud^zf+rn?^ 2 ] . (35) 
J nr. 1 i 



5Y,{s) «5S( S ) 2j Ks 

I use (|17[) for the variation of the integration region, integrate by parts, use the equations of 
motion, to obtain 

= \ ((d,%M<s))f + m^(x{s)f) + J d*x d, (d»%J^) ■ (36) 
The second term can be written as a surface integral. Using (|28|). this becomes 



(5E(a) 2 

- / dE(t)dE(t) <p(t) %d£Gv(t,u) 



(37) 



<5E(s) 



To compute ^^^.s (&) observe that by defining the field <p(x) in the neighborhood 

of the surface by <p(x(s)) = ip(s) and d n ip(x(s)) = 0, the right hand side of ((2*5)) is the flux of 
the vector field 

v»(x) = -<p(x) ^Gx(x,y) , (38) 
Therefore it is possible to use ((T""]) to obtain 

W _ 8 ^ ( „ (s)) * Ge fe(s) ,^ _ / rfE(t) ^ „, (t) JL jg^fatp (39) 



In the first term, I separate the normal and tangential component of the sum over fi and recall 
that <p(x) is constant in the normal direction, and in the second term I use the Hadamard 
formula (|"""""). This gives 

^gf ^ = - V </>(s) V a G s (y( S ), x) - ^( a )D v G? E (y(s),x) 
6 ^ S) (40) 



dZ(t) <p(t) ^G s ( S ,x) S£^G E ( S ,t) . 

But the tangential derivative of the Green function on the boundary surface vanishes, since the 
Green function itself vanishes. So does its Dalambertian, since (— D z + m 2 ) G E (z £ E, y) = 
0). Therefore the first two terms vanish, and using (|28j) again, the result is 

5 %^ } X) = d s n Gx (s, x) d s n %, s (s) , (41) 
7 



where I have used the notation (JT7|) also for After inserting this in (|3TJ) . and using (|HT|) 

the result is 



SE(s) 2 



(42) 

+ / dS(t) <p(t)d a n $(s)d a n d t n Gz(t,s). 



E 



The derivative £]/<fy> (-z) can be easily obtained from the definition . The 

Hamilton-Jacobi equation (|34j) finally reads 



dS{t) d s n $(z) <p(t) d^Gx^s) + ~ (a^( s )) 2 + 

E ^ 

+ ~ [ dE(t)dE(u) ^G E ( S ,t) «G E (s,u) = 



(43) 



E 



2 

Calculating d n $(s) from ((2*5)) it is easy to show that this is an identity. 



5 Generalized Tomonaga-Schwinger equation 

It is possible to write the propagator (|32|) in the form 



1 



W[tp,V] = x/detGs exp (-%,£]]) = exp ^-tr In G E - S^.E] J . (44) 
The variation with respect to the surface (dropping the arguments from the notation) is 

<)U w iJr^(l trlnG x-s) ; (45) 



8E(s) <5E(s) V2 
so that the generalized Tomonaga-Schwinger equation (|24[) becomes 



15(trlnG E ) SS 1 5 2 S 1 £5 6S l f „ ±/ ^ 2 1 2,2/ x 

2 <5S( S ) " 5S(i) = "2^M + 2 M^M^ " 2 (W(S)) ~ 2™ (s) " (46) 

Using the Hamilton-Jacobi equation (|34j) derived in the last section, this equation reduces to 

<J(trlnG s ) 5 2 S 



5S(s) 593(s)5</?(s) 



(47) 



The left hand side of this equation can be calculated again using Hadamard's formula (|33j) 
and J2f 



<5£( S ) " tr 



W7E(zi,* 2 )| t=u= . a , (48) 



££(s) 

while the right hand side can be directly obtained from (|31j) . giving 

1 5^(3)8^8) J ^M*®® V(s) <P® d s n diGx(s,t) = did^G s (z u z 2 )\ t=u=s , (49) 
showing that the equation is satisfied. 



8 



6 Conclusions 



I have derived the generalized Tomonaga-Schwinger equation (J3J) from its functional integral 
definition (j2j), in the case of a free euclidean scalar field. 

An issue that I have not investigated is given by possible restrictions on the shape of the 
region 1Z in using Green functions techniques; for example whether it has to be convex. 

I think that this line of investigation should be relevant in order to develop the program 
presented in [6] to make contact between background independent quantum gravity and 
conventional calculations of particle scattering amplitudes. 

The result derived here undercuts some a priori arguments against the general possibility 
of an equation as the form (J3J), and renders the equation more plausible. However, it is 
of course far from a general proof of the validity of this equation. In the lorentzian case, 
in particular, the distinction between initial and final configurations on the one hand, and 
boundary values of the fields on the other, might appear more substantial. The possibility 
that the variation of the two be governed by the same equation is less evident. The difficulty of 
using functional integral techniques is that integrals that converge in the euclidian oscillate 
in the lorentzian case. If convergence of the integrals can be controlled, I expect that it 
should be possible to prove the validity of the generalized Tomonaga-Schwinger using the 
same strategy as here. 
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Appendix 

A: Hadamard formula 

I sketch here the demonstration of the Hadamard formula (j3.3j) . Consider two surfaces Ei 
and £2, where £2 coincides with Si except for an infinitesimal outward bulge around a point. 
Consider the difference Gs 2 (x,y) — (x,y). This difference satisfies the equation of motion 
If x G S 2 : 



Gs 2 {x,y) - G Sl (x,y)\ 



Gei (x,y) ; 



(50) 



taking y € £ in (|28|) shows that 



d^G s (x,y£Z) = 5W (. 



x 



y) 



(51) 



so that 




(52) 
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The region of S2 where S2 = Si doesn't contribute to this integral because of (|2fjjl. In the 
region where S2 7^ Si, G^ 1 (z,y) is infinitesimal, so it is possible to approximate 

G^{z,y)^d^G^(z,y)8n, (53) 

where n is the normal to Si, positive in the outgoing direction. As a consequence we have 
Hadamard's formula: 

5Gx (x;y) 



<5S(s) 



d^(z,x)d z n G s (z,y) . (54) 



B: The field propagator 



Here I give the explicit expression for the propagator W[<p, S] in the case in which S is 
formed by two parallel planes x 4 = t = and x 4 = t = T. That is, I give the functional field 
propagator of the free euclidean scalar theory. I denote (pi{x) an d <p2(x), the value of the field 
on the initial and final surface, respectively, where x = (x^x^x 3 ) are cartesian coordinates 
on the planes, and I denote the propagator as W[<£>i, (f2, T\. The gaussian integral © can be 
solved by finding the extremal value of the exponent, that is, by solving the classical equation 
with boundary conditions <£(x, 0) = (fi(x) and $(x, T) = if2(x). The solution is 

_ f d 3 kd 3 y i%{x-$) <P2(y) sinhcjt - ^(y) sinhu(t-T) 
1 ' J J (2vr) 3 sinhwT ' 1 ' 

Inserting this in the action and computing the gaussian integral gives 

1 f + °° f +oc f d 3 k 



(56) 



W[ip 1 ,<p2,T\=Nexp{-- / d 3 x / d 3 y / -JL e -*(iM) ^Jp + m 2 
I 1 J -00 J -00 J (2tv) 

f_Min<h&+<h&<hm + coth( ^7^ T) {MylMxl + MylMxl) ) 

V smh(V k 2 + m 2 T) J ) 

The normalization factor N can be found via the Schrodinger equation; the solution is 

(V f d 3 k \ 
N oc exp I — / o lnsinhwT ) (57) 

with V being the spatial volume. For completeness, I give also the propagator in the 
Lorentzian case 



. p+00 f+oo i- j3l, „ r, 
W[i Pl , i p 2 ,T] = Nexp\- / d 3 x / d 3 y / —-3 e -«%-*) \Jk 2 + m 2 

. 1 J -oo J -oo J (27TJ 



sin( v k 2 + m- 2 T) , 



(58) 
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